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Abstract 

We explore cosmology of intersecting braneworlds with induced gravity on the branes. We 
find the cosmological equations that control the evolution of a moving codimension one brane 
and a codimension two brane that sits at the intersection. We study the Friedmann equation 
at the intersection, finding new contributions from the six dimensional bulk. These higher 
dimensional contributions allow us to find new examples of self-accelerating configurations for 
the codimension two brane at the intersection and we discuss their features. 
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1 Introduction 

Brane-world models offer new perspectives for explaining the present day acceleration in purely geo- 
metrical terms, without the need to introduce dark energy [U [21 [3] (for a review see [1]). A celebrated 
example is the Dvali-Gabadadze-Porrati (DGP) model in a 5d spacetime pQ. The brane action includes a 
quantum-induced Einstein-Hilbert action that recovers 4d gravity on small scales. This model realizes a 
so-called self- accelerating solution that features a 4d de Sitter phase even though the 3-brane is completely 
empty. However, so far, only co dimension-one examples of such solutions have been proposed and these 
backgrounds are known to suffer from ghost instabilities [5]. An interesting possibility then is to look 
for other such solutions in higher codimensional set-ups, initially introduced to address the cosmological 
constant problem [6j [Tl [9l [lOl [TT1 [T2] . This might lead to ghost free models [13] (see however [13]). 

In this paper, as a continuation of [15], we consider a codimension two brane that lies at the inter- 
section of two codimension one branes embedded in a six dimensional space. This system was studied in 
the past in the context of standard gravity [16] (cosmological properties were investigated in [E]), and 
Gauss-Bonnet gravity [18] elaborating an idea developed in |20j . The latter was generalized to higher- 
codimensional models in [19J. Models with a generic angle between two intersecting branes were first 
considered in [21] and then further generalised into the so-called Origami- world in [22]. More recently, 
in [15] , we added brane induced gravity terms to this system to analyse the features of a configuration of 
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static branes embedded in a time dependent, maximally symmetric background. We showed the existence 
of new self-accelerating solutions, and of configurations with potentially interesting self-tuning properties. 

In the present paper, we continue the analysis of this system by studying cosmological models, ob- 
tained by the motion of one of the branes through the bulk, in a mirage approach [23] . The energy 
momentum tensor different from pure tension on the branes causes the brane to move and bend in the 
bulk, and induces cosmological evolutions from the point of view of observers sitting on the branes. We 
allow the branes to intersect at an arbitrary angle and to deform in the preferred shape. 

The analysis of gravitational [23] and cosmological [25] aspects of codimension two brane-worlds is a 
subject that is receiving some attention. Cosmology is mainly studied in the context of a mirage approach. 
In higher codimensional brane-worlds, the mirage approach has usually some drawbacks (critically exam- 
ined, for example, in the introduction of [26]), mainly due to fine-tuning relations that the brane energy 
momentum tensors must satisfy. These are usually associated with the fact that an analogue of Birkhoff 
theorem does not hold in this case, in contrast to the codimension one case. In codimension one case, 
this theorem ensures that a system composed by a homogeneous and isotropic brane, moving through a 
static higher dimensional space, fully catches all the relevant time dependence of the system [27] . In our 
case, this is not true: generically, a moving higher co-dimensional brane induces time-dependent effects 
in the bulk \18\ [28] . In order to avoid the time dependence in the bulk, one must impose a static ansatz 
for the bulk geometry, and this is reflected on fine-tuning relations between matter on the brane and 
in the bulk. Nevertheless, it remains the most direct approach to study cosmological aspects of these 
models analytically. 

The most interesting problem in this system is the isolation of the six dimensional effects in the induced 
Friedmann equation on the codimension two intersection. As we will see, the Friedmann equation at the 
intersection receives contributions due to induced gravity terms on it which ensure the recovery of normal 
4d cosmology in the relevant regimes. Moreover, there are terms coming from induced gravity on the 
codimension one branes, of the typical DGP form [2]. Finally, and more interestingly in our framework, 
the Friedmann equation contains also contributions that come from the six dimensional bulk. They vanish 
in the limit in which the branes intersect at a right angle, but for generic brane configurations they play 
an important role for the cosmological evolutions. Indeed, they can provide the late time acceleration, 
regardless of the energy content of the codimension two brane, generalising the self-accelerating branch 
of the codimension one DGP model to higher codimensions [2]. This fact has been realized already in 
[15] . but the present analysis is more general because we do not impose the maximal symmetry on the 
branes under consideration. By properly choosing the embedding for the codimension one branes, the six 
dimensional effects at the intersection can depend on the inverse of the induced Hubble parameter, and we 
will analyse the consequences of this in our discussion. Another peculiar feature of our construction is that 
six dimensional contributions to the Friedmann equation are also associated to the non-conservation of 
the energy density at the intersection. During the cosmological evolution, the energy density indeed flows 
from the codimension one to the codimension two branes, unless the codimension one branes intersect 
with a right angle. 

This paper is organised as follow. In Section [21 we will present the general formalism that is necessary 
to study cosmological properties of the systems we are interested in. In Section [3j we apply this formalism 
to a particular embedding for the codimension one branes, and in Section [J] we study in some detail 
cosmological solutions derived from this embedding. Then, in Section (5J we study applications of these 
cosmological solutions to some interesting situations. We conclude in Section [6] 



3 



2 The general formalism 



2.1 The model 

We consider a system of two intersecting codimension one branes embedded in a six dimensional space- 
time. They intersect on a four dimensional codimension two brane, where observers like us can be 
localised. We take an Einstein-Hilbert action for gravity in the bulk and we allow for induced gravity 
terms on the codimension one branes, as well as on the intersection. Besides gravity, we allow for a 
cosmological constant term in the bulk, A#, and for additional fields localised on the branes described 
by general Lagrangians L's. The general action takes the form 



+ £ d^xV^ (¥±Rn + Ln\ (1) 

where S n = Hi denotes a three- brane at the intersection between all codimension-one branes £j. We 
can have different fundamental scales in the different regions of the space, Mq, M§ t j, and M4. The induced 
gravity terms could be generated, as it was proposed in the original model, by quantum corrections from 
matter loops on the brane. It is also interesting to note that induced curvature terms appear quite 
generically in junction conditions of higher codimension branes when considering natural generalisations 
of Einstein gravity |29[ [30] as well as in string theory compactifications [31], orientifold models and 
intersecting D-brane models [32] . 

The six dimensional bulk is characterised by a maximally symmetric geometry 



ds 
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A 2 (t, z l ,z 2 ) [rj^dx^dx" + b kh dz k dz y 



= 1 T fFTT I i ■ ( 2 ) 
1 + Hi + kiZ 1 



The parameters H and ki appearing in the warp factor A satisfy the following relation 

H 2 — k\ — fef, (3) 
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in order to solve the Einstein equations in the bulk. 

We embed a moving and a static codimension one branes (£2 and Si respectively) on the background 
given by ([2]). The moving brane £2 is characterised by an embedding 

= (i, x 3 , Zxiun), Z 2 {t,uo 1 )) . (4) 

Here, W\ is an embedding coordinate. In the following, for simplicity, we will demand that the 
intersection with the other brane lies at the position w\ = 0, and that the function Z\ does not depend 
on time. The vectors V tangent to £2 are given by (we introduce indices on the left of V: they indicate 
which brane we are talking about) 

FlY M 

(2) K)=^ * a = (t, X \ Wl ), (5) 
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where is (proportional to) the velocity vector 

= (1, 0\ 0, Z 2 ) = X M . (6) 

The other four vectors are 

m V M = (0,^,0,0), (?) 

(2) < l} = (0, 0\ Z[, Z> 2 ) = X' M . (8) 

The normal vector to the brane is thus given by the conditions 

n M V$ = 0, Vo. (9) 

Orthogonality with respect to the i vectors simply removes from um all its 3-dimensional space- like 
components. Imposing orthogonality w.r.t. ^'VM\ we then find 

with 



N = dZ? + Z»-Z*Z». (11) 



X% = (t, x 3 , 0, z 2 ) , (12) 



Doing exactly the same steps for the static brane Ei, with embedding 

rM 

the vectors tangent to the brane, ^V^, are immediate to find. And the normal is simply 

n { l) = A (0, 3 , 1, 0) . (13) 
2.2 An useful change of coordinates 



Proceeding identically as in the static case |15[ [22] , it is useful to change a frame in order to impose the Z 2 
symmetry in the case of a general angle. We go to coordinates parallel to the branes (z±, z 2 ) — ► (z±, z 2 ), 
where 

dz k = n (fe ) • dz . 
One obtains two two-vectors lm and 1(2) parallel to the branes: 

!(D = ^rC 2 !^)- 1(2) = ^GU), (14) 

and dz = dz k . Then the components of the vectors parallel to the moving brane become 

BX M 

m ^ = ftf " *■ z; ' °) • (15) 
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and 



(0) 



1, 3 , : 



(16) 



Notice that the consistency relation 



d 2 X§ _ d 2 X^ 



in our case implies the condition 

dwi \ M 

The normal to the moving brane becomes in these coordinates 




(17) 



n 



(2) 
M 



AT 



Z2 Z[ 

'~at 



, o 3 , 0, 1 



(18) 



where T is a normalization factor that we will fix once we define the six dimensional metric. 

In order to proceed, we must take into account that the branes are fixed points of Z2 symmetries. We 
focus our analysis on the moving brane £2 in order to compute Israel junction conditions at its position: 
the analysis for the static brane £1 can be performed along similar lines. The Z2 symmetry acting on the 
static brane £1 implies the invariance of the 6d metric under z 1 — > —z 1 , that can be obtained replacing 
z 1 - Iz 1 !. 

After the change of frame, imposing the Z2 symmetry, the six dimensional metric becomes 



7,i 



1 

Z* 



with inverse 
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Z? + z'i 

M Z' 2 sign(zi) 
N 2 

—J\fZ' 2 sign(zi) 



J\f Z' 2 sign(zi) 
—NZ' 2 sign(zi) 



(19) 



Z' 2 + Z' 2 



where 



1 + 



?' 2 
Z' 2 



(l - sign 2 (zij) . 



(20) 



(21) 



With this information, we determine the normalization factor T by requiring that the normal 
has a unit length. Then T is determined as 



T 



CM 
~2~' 



with J\f defined in formula (jlip . C in (|2 lfl . and 

Q 2 



Z? + Z' 2 2 



c 2 z' 2 z 2 2 . 



(22) 



(23) 
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On the other hand the induced metric is invariant under bulk reparametrisation, and thus reads 
<&§,Ea = A 2 (*> w l) ~ (l - ^l) ^ 2 + + (^f + Z 2 ) dwj 

+ 2Z 2 Z' 2 sign(z 1 )dtdw 1 . (24) 

The inverse induced metric is given by 
with 

V Z 2 Z' 2 signih) 1-2% J 
and Q defined in eq. (123p . At the intersection, the four dimensional metric is given by 

ds 4 = A 2 (t)[-{l-Z 2 2 )dt 2 + dx 2 }. (27) 

2.3 Extrinsic curvature 

Given all this information, one can compute the components of the extrinsic curvature at the position of 
the brane T, 2 , using the general formula 

K mn = V^V^Vmun . (28) 

Since the expression for K mn is invariant under bulk reparametrisation, to evaluate the right hand side 
of the previous expression one can use the six dimensional metric in the original frame, or in the frame 
parallel to the brane. 

The calculation of the regular part of K mn is easier to work out in the original frame. The non 
vanishing components are the following 

K 00 = ~4-A Z 'i " T7 f 1 " Zi) *) > (29) 



* { z»z>-z»zi l)+ * 

A _. • , . A 2 



K m = - (Z'(Z' 2 - Z»Z[) + (Z^ 2 + Z 2 ' 2 ) IC(w, t) , (30) 



K 0wi = -sign(z!)jj Z[Z' 2 + sign(zi) -jj K.(w,t) Z 2 Z' 2 , (31) 



A 2 5jj %i) 
j = Jf ' 



where 

K{w, t) = hZ' 2 - k 2 Z[ - HZ 2 Z[ (33) 

In addition, the component K Wim of the extrinsic curvature may contain terms localised at the inter- 
section due to the presence of the sign functions in the six dimensional metric. Let us then look for the 
singular pieces of the extrinsic curvature 

K ab \ sing sing ■ 

(34) 
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This quantity is much easier to calculate in the tilted reference frame. There are, a priori, two classes 
of contributions to the singular pieces. The first one comes from partial derivatives acting on njy, due 
to the sign function included in T. However such a contribution is proportional to njy itself as the sign 
function only appears in the prefactor of njv/i and thus vanishes due to the orthogonality between um 
and V, A '{. 

(a) 

Then 

K ab \sing — V( a )V(b) n R ^ MN I sing-, (35) 

and since has no singular part, one is left with T^ N . Its only singular component is r? 2 ^ = g 22 d\g\2 
since they involve derivatives of the sign function. Then, in the end, we find that 

K WlWl \ sing = -AK, 1 2 Z' 2 —sign(z\) , (36) 



C 2 M 2 \dz! 

is the only singular component of the extrinsic curvature. Notice that the previous expression contains 
products of distributions, since the C contains squares of sign functions. In order to dispel any doubts 
about how to define such singular expression, it is convenient not to set - st g^ Zl - = 2 - , but instead 
maintain the derivative of the sign function. Later, in the specific examples we will discuss, in order 
to extract the singular terms localized at the intersection we will perform an explicit integration on a 
small interval around the singularity. The result of this integration will provide the value of the various 
quantities localized at the intersection. 

For the static brane Ei, the extrinsic curvature is simply given by 



Kt = -hS a b . (37) 



2.4 Junction conditions 



The previous expressions for the extrinsic curvature are important in order to obtain the equations that 
govern the induced cosmology on the brane. They are dictated by the Israel junction conditions 



K, 



ah 



2 [K ab - Kh ab ] = -j^(S ab + S l a f), (38) 



where [X] = (X(S2,+) — X(£2,-))/2, while the induced codimension one brane metric is h a \,. The 
extrinsic curvature tensor evaluated on £2 is given by K ab = h Al 'h^V ' m^n with K = K a a , and energy 
momentum tensors relative to matter localised on £2 1 appearing on the right hand side of (|38p . are 
calculated in the usual way: 

bab -, — 1 , (39) 

V _ft (2) dtl (2) 

S l : b c = -<5(£!)« 



y/~ h (2) Sh n 



'-^«^- (40) 

In our model the localised energy-momentum tensor also includes contributions from the induced gravity 
terms. The last quantity S l ° b c denotes energy momentum tensor that is localised on the intersection £ n 
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between the branes. Notice the presence of the factor y/in/^(2) that renders the expression covariant 
with respect to the metric at the intersection. 

In the previous discussion, we learned that the only singular term of the extrinsic curvature for 
the brane £2, that is localised at the intersection, is contained in K W1W1 . This implies that the six 
dimensional contributions to the energy momentum tensor must be proportional to the induced metric, 
Suv = f( xlJl ) hr\/j,vi f° r some function /. We still do not know whether this function / is a constant (in 
which case, it corresponds to a pure tension) or not, since we do not know whether the conservation of 
the energy holds at the intersection or not. The Codazzi equation holds in this case0 

V a K% = =► V a S% = 0, (41) 

which means that there is no exchange of energy between the bulk and the codimension one branes. 
But the previous relations may contain singular terms, associated with an exchange of energy between 
the codimension one and the codimension two branes. This is indeed what generically happens, and we 
will encounter an example of this phenomenon in Section O Singular terms in the first of the previous 
formulae can appear if has singular pieces, or if some of its components become singular when covariant 
derivatives act on them. This possibility occurs when the angle between the brane is not right: then the 
component -Kg"" 1 IS non-vanishing at the intersection, and, being proportional to the sign function (see 
eq. (|3ip ). it normally generates an additional singular term. 



3 Applications 

We then consider a system with the static brane Si, and the moving brane £2. Here, £2 is free to move 
and bend arbitrarily. We assume that the induced gravity term on the moving codimension one brane 
vanishes: Afg 2 = for simplicity. For this brane, we take an embedding 

X M = (t, x 3 ,Z 1 (w 1 ),Z 2 (t,w 1 )), (42) 

with 

Z\ = w\ cos a(wi), (43) 

Z 2 = Z2(t,wi) + wi sina(wi), (44) 

where we wrote the two functions Zi in terms of the auxiliary functions z 2 and a. We demand that 
these functions are continuous with respect to the variable Wi, at the position of the intersection w\ = 0, 
and to avoid subtleties related with the reflection symmetry at the intersection we require smoothness 
conditions z' 2 (t,0) = 2%(t,0) = a'(0) = 00. From these definitions, we have 

Zi =0, Zi = cos a — w\a' sin a, 

• (45) 
Z 2 = z 2 , Z' 2 = z' 2 + sin a + w\ct cos a. 

Notice that the previous embedding satisfies the constraint (JTTj) ; for a time-dependent angle (I17j) would 

1 The RHS of this formula vanishes because both the bulk and the static brane have maximal symmetry [33] . 

2 Asking that only the first derivative vanish at the intersection may be enough to ensure sufficient smoothness to render 
the system well behaved. In fact later we will briefly mention a situation where a non- vanishing second derivative z'^it, 0) ^ 
can turn out to be useful. 
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instead imply a(t, W\ = 0) = 0. 

The junction condition on the static co-dimension one brane gives the tension Ai as 

Ai = GMl^H 2 - fef ) + 8Af 6 4 A;i . (46) 

The induced metric on the brane £2 is obtained by plugging the previous expressions in (|24p . The 
complete calculation of the cosmological behaviour for the moving codimension one brane is complicated 
as the brane is inhomogeneous, but it can be obtained straightforwardly from the general formulae ([291) 
(|32p. In the next subsection we discuss some of its properties that are useful when comparing them with 
cosmology on the codimension two brane. 



3.1 Cosmology on the moving four brane 



The cosmological evolution on the moving brane £2 is complicated by the fact that its induced scale factor 
and energy momentum tensor must be inhomogeneous, in order to satisfy Israel junction conditions (|29l) 
(|32p : such conditions require some off-diagonal components of the energy momentum tensor S a i, to be 
non vanishing. In what follows we will only need the explicit form of the codimension-one equations 
evaluated at the intersection. We thus concentrate on such a limit where the form of the needed energy 
momentum tensor is the following: 



( -p 2 \ 

p 2 

p 2 

p 2 

V x p 2 ) 



at w\ = . 



(47) 



Then the junction conditions impose the following relations 

p 2 = -SMlN- x 1C, 
p 2 = 8Af 6 4 (AA~ 1 /C + -N^A^hcosa) 
X = 2Mq A~ x Z2Z2 sin a cos a, 



(48) 
(49) 
(50) 



where we define N = \J\ — cos 2 az\ and K, = k\ sin a — (&2 + HZ2) cos a. 



3.2 Cosmology at the intersection 

We start from discussing the contributions from the brane £2 to the codimension two brane. At the 
intersection, characterised by w\ = 0, the induced metric is straightforwardly extracted from the five 
dimensional one and is simply given by 

ds\ = A 2 (0,t) I - [1 - if] dt 2 + dx*l} = -dr 2 + a(r) 2 dxl, dr = A(0,t)yjl- z%dt. (51) 

Then the induced Hubble parameter is 

1 da(r) H + k 2 z 2 , . 

a(r) dr _ z | 

In order to find the Friedmann equation at the intersection, we have to extract the singular part of 
the Israel junction conditions for the codimension one branes. This singular part receives contributions 
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from the energy momentum tensor localised on the codimension two brane (containing also the induced 
gravity terms at the intersection), from the induced gravity terms on the codimension one branes, and 
from singular contributions of the extrinsic curvature terms. The final contribution represents the most 
interesting feature of our model since it corresponds to six dimensional contributions to four dimensional 
physics. We start our discussion with their evaluation. 

Recall that only singular part on the extrinsic curvature for brane £2 is contained in the (w\,Wi) 
component. From its expression one straightforwardly obtains 

^ ( "" 9> = ZvW (1 " il) ( Jr 9n(2l) ) «• <53) 

as the only singular piece, and it is easy to see that its contribution to the energy momentum tensor 
at the intersection results proportional to the induced metric (recall the definitions of T, jV, C and Q 
respectively in eqs. (]22|) . (fTTI) . (f2"T1) and (123]) ). Let us, for example, calculate its contribution to the energy 
density at the intersection. Plugging the previous expressions inside the Israel equations, one finds the 
relation 

5p ( d . \ T tana r — / d . \ 

ami \M sl9n{zi) ) = -cWq Vl ~ z [w^ 9 ^ ) ' (54) 

where 5p indicates the contribution to the brane energy density, due to purely six dimensional effects. 

The previous equation contains sign 2 (z\) functions inside the expressions for Af, Q, and C. The safest 
way to handle them is to integrate both sides of (|54l) along an infinitesimally small interval centered at 
the origin, furnishing the value of the contribution to energy density localized at that point. Performing 
the integration we get the expression 

5p f ~ f +1 dsign(zi) 



- tana-wl — zz / — —— ^ — . (55) 

2M 6 4 V 2 J_ x C Q 2 y ' 

After a few calculations, the previous expression can be written as 

5p Sina f+1 dx 1 1 (56) 



2M 6 4 V 1 - M, J- 1 \/l - sin" a a? 1 + ^.." x 

By performing the integration, we get 



5p = 4Afi arctan 



tan a 



0" 



(57) 



Notice that, in the limit of static brane Z2 = 0, this contribution is proportional to the angle a, exactly 
as happens in the case of codimension two conical singularities. 

Proceeding with our calculation, we can determine the contributions to the intersection from the 
induced gravity terms on Si (recall that we have chosen M5 2 = so there are no induced gravity terms 
on the moving brane £2). We find 

joying = _6(fe2+£j). (58) 

y/l -Z$ 
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Putting all this information together, we find the following equation relating energy density to geo- 
metrical quantities (here p indicates the total energy density at the intersection) 



99 ? ko + Hzo A ( tan a \ , „ 

p = 3M%H 2 + 6M| 2 + 4M 6 4 arctan - , (59) 

Vi -i 2 2 Vv 1 -^/ 

that can be interpreted as the Friedmann equation for an observer localized at the intersection. The 
induced gravity terms on the codimension one branes do not induce a violation of the continuity equation 
at the intersection because the intersection can be seen as codimension one object from the point of 
view of the four branes and then the properties of the Israel formalism for the codimension one brane 
ensure the conservation of energy (see eq. (|4ip ). On the other hand, the last, six dimensional term in 
Eq. (|59j) is explicitly time dependent, while we know that it appears as a tension term in the effective 
energy momentum tensor at the intersection. This is because it is proportional to the induced metric. 
This indicates that this term is likely to be associated with a violation of the continuity equation at the 
intersection. 

This issue can be understood by re-considering the Codazzi equation: 

V„Kjf = (v A/ ^) (rCS) + (vM4 f ) (Sing) = 0, (60) 

where V is the covariant derivative with respect to the five dimensional metric on the codimension 
one brane. From the previous formula, we learn that both the regular and singular parts must vanish 
simultaneously. However, it can happen that the covariant derivative induces singular contributions when 
it is applied to certain components of K^f by taking derivatives of sign functions. This is indeed what 
happens in our case. Consider the case N = 0. The singular part of the previous formula tells us that 

±(K%)" n ' = (V M K^) sin9 , (61) 

where the left hand side contains the singular term associated with K WlWl \ s i ng , while in the right hand side 
the singular terms come from the covariant derivatives. But the piece in the left hand side corresponds 
precisely to the term associated with the six dimensional contribution at the intersection. Thus the 
six dimensional contribution to the Friedmann equation on the intersection does not satisfy the energy 
conservation and there is an exchange of energy from codimension two brane to the higher dimensional 
space. 

In the light of this fact, one expects that the conservation of energy at the intersection does not hold. 
Instead, one finds the continuity equation 

A d I tana \ 

p + 3H (p + p) = 4M 6 4 — arctan . (62) 

dr \y/\-zl) 

Hence, the conservation of energy is ensured only when a vanishes, or when z% is constant. 

We close this section by summarising the equations that govern cosmology of the codimension one 
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branes and the intersection; 



^ = H^-kl-kl (63) 
Ai = 6M'i(H 2 -A;|) + 8M 6 4 fci, (64) 
P2 = —8Mq — — — (ki sin a — (k 2 + Hz 2 ) cos a) , (65) 



cos^ az. 



on ,4 2 2 z 2 smacosa 
X = 2Me A(0,t)(l- cos a%)W ^ 

p = 3M|^ 2 + 6M| 1 ^±£f + 4M 6 4 arctanf^=Y (67) 

where A(0,t) = 1/(1 + ift + ^^(i)), A# is the bulk cosmological constant, p2 is the energy density on 
the moving codimension one brane at wi = 0, x is ( w i > t)-component of energy momentum tensor on the 
moving codimension one brane at w\ = and p is the energy density at the intersection. The Hubble 
parameter at the intersection is given by Eq. (|52p . The energy conservation at the intersection is given 
by Eq. (f62l) . In the following we put Mi 5 = M 5 . 



4 Cosmological solutions 

In this section, we discuss the property of the cosmological solutions by focusing on the Friedmann 
equations on the moving four brane and at the intersection. 

4.1 The branes at a right angle 

We first consider the simplest case in which the branes are at a right angle a = and H = 0. In this 
case, the energy momentum tensor on the moving four brane becomes the perfect fluid and there is no 
energy flow x = 0. Using the cosmic time r, the 5D metric is given by 

dsl = -dT 2 + A\w llT ) (dxl + dwl), A(w 1 ,t) = -— — — . (68) 

Although the scale factor depends on wi, the Hubble parameter in terms of the cosmic time is independent 
of w\ and given by Eq. (|52p . By expressing Z2 in terms of the Hubble parameter H, we get 



p 2 = 8M^H 2 + kj. (69) 
Since there is no energy flow, the energy density is conserved 

d T p2 + 4H(p 2 +p 2 ) = 0. (70) 
At the intersection, the Friedmann equation is given by 



/ H 2 

p = 3MjH 2 + 6Mik 2 Jl + -p-, (71) 

and the standard continuity equation holds since tana = 0. Notice that the static brane gives a contri- 
bution of the 5D DGP form. 
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Since the Hubble parameters are equal in both the equations (|69p and (|7ip . by expressing H as 
a function of the energies in the two cases and equalling the results, they will imply a fine tuning 
relation between the two homogeneous energy densities p2 and p. Geometrically, this is because when 
the codimension one brane S2 moves through the static bulk, it completely controls the dynamics of the 
brane S n that sits at the intersection with Si. Then, S n can only follow the motion of S2, without an 
independent dynamics on its own. The problem becomes clearer by the fact that the energy density and 
the Hubble parameter on the moving codimension brane do not depend on the coordinate w±. Then, the 
energy density at the intersection actually fixes all the properties of the energy density on the moving 
brane S2, including its equation of state. 

There is a simple way out of a part of this problem. The fine-tuning we found is so strong because 
we demand that the moving brane S2 keeps a straight shape - that is, it cannot deform along the z\ 
direction. Suppose however that we allow the moving codimension one brane to be free to deform and 
bend, forming a non-trivial angle a with Si which explicitly depends both on z\ and t. Then, the 
energy density on S2 will explicitly depend on z\. This implies that, although the energy density at the 
intersection must equal the energy density on S2 calculated at Z\ = 0, nevertheless this fine-tuning is 
ameliorated with respect to the previous case. Indeed, it involves only the quantities calculated at the 
intersection. 



4.2 Arbitrary angle between the branes 



Now let us consider the case in which In this case, the Hubble parameter at w\ = is given by 

Eq. (|52p . Again taking H = for simplicity, the Friedmann equation on the moving four brane at the 
position of the intersection is given by 



P2 = 8Mg (fo cos a — ki sin a) 



H 2 + k 2 



H 2 sin 2 a + k\ ' 

On the other hand, the Friedmann equation at the intersection is given by 



(72) 



3M%H 2 + 6Mf k 2 



I ff2 

1 + — T + 4M@ arctan 



tana* / 1 + 



k\ 



(73) 



When the angle a vanishes and k 2 remains finite, we recover the results of the previous subsection. 
Notice also that comparison between Eqs. (172p and (1731) imposes a fine-tuning relation between p and the 
energy density pi of the codimension one brane, when evaluated at the intersection. Nevertheless, this 
fine-tuning is much milder than the one we met in the previous subsection. This fine-tuning is associated 
with the restrictive Ansatz we have chosen for the bulk metric. 

The form of the previous Friedmann equation is quite complicated to study with full generality. While 
the second term on the right hand side of Eq. (I73p corresponds to the well-known DGP-like term, the 
last term in the right hand side of Eq. (|73p is less standard, and is associated with six dimensional 
contributions. 

Notice that, in the limit of large H, this term approaches the constant value 4Mg sign(k2), and 
then can help to drive acceleration in this regime. We will discuss in the last part of the paper more 
general situations, where six dimensional contributions can provide sources of acceleration for the induced 
cosmology at the intersection. 
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The continuity equation is given by 



p + 3H (p + p) = 4Mi — arctan 

or 

so, under the assumption that hi is non zero, the conservation of energy is ensured only when a vanishes, 
or when H is constant. For a^O, it is also necessary to have the energy flow on the moving codimension 
one brane, x, which also breaks the conservation of energy on the moving codimension one brane. Then 
we can understand that the energy at the intersection is transmitted to the moving codimension one 
brane. 



tanau / 1 + 



IP 

^2 



(74) 



5 Applications 

In this section, we derive some interesting consequences of the cosmological equations we discussed in 
the previous section. In the first two subsections, we examine some of the cosmological properties of the 
solutions we discussed in [15] in the present context. In the last subsection, we will instead derive a new 
selfaccelerating configuration in which the codimension one branes are not maximally symmetric, a case 
that we did not discuss in our previous work. 

5.1 Self-tuning solutions 

Here we re-examine the selftuning solution presented in [15]. We take k\ = 0, ki = 0, M5 \ = and 
i = 0. Then we have 

A 4 = 3M%H 2 + 4M 6 4 a, H 2 = (75) 

where H = H. Then the expansion rate does not depend on A4. The self-tuning mechanism consists on the 
fact that if we change tension A4, a changes so that induced cosmology remains the same. Unfortunately, 
our embedding is not well suited to study the self-tuning property of the solution as a = constant is 
actually imposed by hand. A possible way out would be to consider a situation in which z"(t, 0) 7^ 
at the intersection. This would generalise (|17l) with new pieces that would not necessarily impose that 
a(t, 0) = 0. It would be nice to study in more detail this kind of generalisation to understand whether 
it can be compatible with the reflection symmetries of our system or not. It is important to understand 
whether the eventual self-tuning property would be compatible with the recovery of small scale 4d general 
relativity on the intersection, in order not to contradict big bang nucleosyntesis and other cosmological 
tests. In order for this last tricky issue to be solved, it seems to be necessary that the dynamical angle 
reacts only to the vacuum energy density component of the localised matter on the intersection: a priori 
this is rather counterintuitive. However, a few observations are in order here: it is well known that 6d 
brane worlds with conical singularities treat tension-type of matter on a completely different footing with 
respect to a generic fluid (uj 7^ —1) [34j . In fact, also in the present setup the six-dimensional contribution 
to the 4d stress tensor has a tension-like structure; moreover as we showed, a generic fluid localised on 
the intersection does not seem to render the bulk geometry singular as opposed to what happens in the 
thin conical setups [33]) but violates the conservation of energy. It is therefore not excluded that the 
self-tuning might be at work here. 
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5.2 Self-accelerating solutions with H ^ 



Here we re-consider the self-accelerating solution presented in [T5], for maximally symmetric configura- 
tions. We assume there is no cosmological constant nor matter in the system Ab = p = p 2 = Xi = with 
z = 0. Then we get 










H ki k 2 , 
3Mf # 2 + 6Mlk 2 + 4M|a, 
k± sin a — k 2 cos a, 
6Ml(H 2 - kf) + 8M$ki. 



(76) 
(77) 
(78) 
(79) 



If k\ < and k 2 < 0, there are non-trivial solutions for k±, k 2 , ce and H. The solution is roughly given by 

R ~ld A ' M B 3 ~M 4 M|. (80) 

in accordance with what we found in the previous paper. Once a, k 2 and H are fixed, Eqs. (|62p and (|67p 
determine the cosmological dynamics with p without ambiguity. The resulting cosmology is complicated 
due to the non-conservation of energy. Instead of dealing with this complicated case, we will discuss a 
simpler, different situation in the next subsection. 



5.3 Self-accelerating solution with H = 

We consider the case H = 0, k\ = 0, and k 2 = —tiji \ sina| Mq, for some positive constant (3 and e 2 = ±1. 
To make a more direct comparison to analog studies in the codimension one case, we focus on a regime 
in which the quantity H 2 satisfies the condition 



H 2 > /3 2 sin 2 aMi, 



(81) 



and we will later check in which cases this relation can be satisfied in our context. The continuity equation 
becomes 



p + 3H(p + p) 



d 

-4 ei Ma — — arctan 
or 



H 



j3 Mq cos a 



(82) 



with ei = — i 



Jr = ±1. Notice that ei = 1 means that we are taking negative values for the angle a: 
in the conical case, it would correspond to take an excess angle for the conical singularity. 
On the other hand the Priedmann equation acquires the form 



P 



3 Mj H 2 - 6M| e 2 H - 4 e x M| arctan 



H 



(3 Mq cos a 



(83) 



with e 2 = ±1 corresponding to the usual DGP choice of branches for the codimension one brane. The 
first term in the right hand side is dominant at large H, and ensures the correct four dimensional form 
for early time cosmology. The second term is the typical DGP contribution, while the third term, a 
six dimensional effect, is less standard as we discussed before. The previous expression can be easily 
rewritten as 

Mf\ 2 Ml 4eiM 6 4 H 
H -e 2 -rS; I - t4- — ,o arctan 



3M 2 



Ml 



M\ 



3M 2 



(3 Mq cos a 



(84) 
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from which we obtain 



H 




Ml 



+ ei 



i_M£ 
3M| 



arctan 



H 



(3 Mq cos a 



(85) 



by imposing that the quantity inside the square root is positive. Now, the choice e\ = £2 = 1 corresponds 
to the standard DGP self-accelerating branch, and the quantity inside the square root is always positive. 
It implies that, even when p vanishes, the Hubble parameter satisfies the inequality 



H > 



Ml 
Ml 



1 + 4/1 + 



4 M| M? 
3M| 



arctan 



II 



(3 Mq cos a 



(86) 



and so we find a lower bound for H, as in the well-known self-accelerating branch of DGP model. 

This case is very similar to the standard five dimensional case, since the acceleration is mainly driven 
by the effects of the codimension one brane. It is however also possible to study the case in which M5 = 0, 
to understand whether six dimensional effects provide acceleration by themselves. We will focus on this 
case, choosing t\ = 1, for the remaining discussion. Then, the continuity equation f)82f) can be formally 
integrated as 



P 



Po 



a 



-3(l+w) 



exp 



X(H) 
X(oo) 



dr 



3A/, 
AM, 



^(3 2 cos 2 a tan 2 x 



(87) 



defining 



X(H) = arctan 



H 



(3 cos a Mq 

and introducing the constants po, ao (corresponding to the quantities evaluated at a fiducial time) and 
calling w = p/p. This solution shows that, in the limit H — > 00, one recovers the usual relation between 
energy density and scale factor because in this limit the last factor in the previous equation becomes 1. 
Using the Friedmann equation, it is straightforward to get the following relation for the acceleration: 



a 
a 



H + H 



1 



2M| 



1 



+ w 



8M 6 4 



arctan [H/ ((3 cos a Mq)] 



3M| H 2 



4 M| arctan [H/ {(3 cos a M 6 )] 



so we have the acceleration when 



u < 



+ 



8Mi 



X(H) 



3M 2 H 2 



A(3M*X(H) 



(89) 



(90) 



Then for small H we learn that the six dimensional contributions help to provide the acceleration and it 
can be achieved even when lj > — 4. 

To conclude, we discuss the late time cosmological evolution for our system. In an expanding universe, 
at late times p — > 0; then H approaches a constant value, given by the solution of the equation 



4M| 
3M| 



arctan 



H 



(3 cos a Mq 



0. 



(91) 



Let us study two limiting cases, compatible with the condition (|8ip . in which equation (|9ip admits simple 
solutions. The first is the case in which a is not too small, from which applying (|8ip in (|9ip we get the 
equation 

r2 



4M 6 4 



3M| 



(92) 
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In this case, H can be rendered small by choosing M4 sufficiently larger than Mq. On the other hand, 
condition (|8ip requires Mq (3 M4, so the constant j3 must be chosen correspondingly small. Notice 
that relation (|92|) is similar to the one we already met in eq. (|8U|) . although in the present case we have 
M 5 = 0. 

The second possibility is to consider a extremely small, and at the same time focus on a regime in 
which H -C (3Mq. The the solution of (|9ip provides a different relation between H and the mass scales: 

4M| 

H -jmi' (93) 

and also in this case we can get a sufficiently small value of H by choosing M4 much larger than Mg. 

Notice that the present examples of self-acceleration require a different analysis, with respect to the 
ones discussed in the previous subsection originally found in [15]. This is because, for the particular choice 
of our embedding, the codimension one branes do not need to be maximally symmetric, nor empty. It 
would be nice to understand whether in this case ghosts are present in the low energy spectrum, and if 
so how do they manifest themselves. 



6 Conclusions and Open Issues 

In this paper, we explored the cosmological features of a codimension two brane with induced gravity 
terms, sitting at the intersection between two codimension one branes in six dimensions. We found that 
the cosmological expansion at the intersection is controlled by contributions coming from the codimension 
one branes, and from the six dimensional bulk. We first showed that the effect of the codimension one 
branes on the Friedmann equation at the intersection have the well-known DGP form. Then, we learned 
that six dimensional contributions are much less standard. They can have an important role for late time 
cosmology providing a new source of the geometrical acceleration, controlled by the angle between the 
branes. At the same time, they are also associated with a violation of the energy conservation at the 
intersection, allowing a flow of energy between the codimension two brane and the higher dimensional 
space. We discussed consistency relations that matter on the codimension two brane must satisfy and 
the connection with the choice of energy momentum tensor localised on the codimension one branes. 

The main aim of this work was to formulate a general and powerful formalism based on the approach 
of mirage cosmology that can be used to study cosmological solutions in this and similar models, and 
to apply it to a couple of representative examples. Due to the fact that the codimension one branes 
that intersect with general angle are not homogeneous and isotropic, a numerical analysis is likely to be 
needed in order to analyse in full details the cosmological evolution of this class of models. As a natural 
continuation of the present work, it would be interesting to study the low energy effective action for 
the light modes associated with our brane configurations. This analysis would be necessary in order to 
investigate whether ghosts are present in the spectrum of the low energy theory, and, if so, whether they 
can be eliminated with a mechanism similar to the one of [13] . We leave these issues to future work. 
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